Single wall carbon nanotubes as coherent plasmon generators 
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FIG. 1: (Color online) The geometry of the problem. 
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The possibility of low-energy surface plasmon amplification by optically excited excitons in small- 
diameter single wall carbon nanotubes is theoretically demonstrated. The nonradiative exciton- 
plasmon energy transfer causes the buildup of the macroscopic population numbers of coherent 
localized surface plasmons associated with high-intensity coherent local fields formed at nanoscale 
throughout the nanotube surface. These strong local fields can be used in a variety of new optoelec- 
tronic applications of carbon nanotubes, including near-field nonlinear-optical probing and sensing, 
optical switching, enhanced electromagnetic absorption, and materials nanoscale modification. 

PACS numbers: 78.40.Ri, 73.22.-f, 73.63.Fg, 78.67.Ch 



I. INTRODUCTION 

Single wall carbon nanotubes (CNs) — graphene sheets 
roUed-up into cylinders of ^1 — 10 nm in diameter and 
~ 1 /im up to ^ 1 cm in length [H-Q — are shown to be 
very useful as miniaturized electromechanical and chemi- 
cal devices Q, scanning probe devices_pi,ii|, and nanoma- 
terials for macroscopic composites [3-01 • The area of their 
potential applications was recently expanded to nanopho- 
tonics [10-15| after the demonstration of controllable 
single-atom incapsulation into single- walled CNs [l6j . 
and even to quantum cryptography since the experimen- 
tal evidence was reported for quantum correlations in the 
photoluminescence spectra of individual nanotubes [iTj . 

The true potential of CN-based optoelectronic device 
applications lies in the ability to tune their properties in a 
precisely controllable way. In particular, optical proper- 
ties of semiconducting CNs originate from excitons, and 
may be tuned by either electrostatic doping [H, [3 ' 
or via the quantum confined Stark effect (QCSE) by 
means of an electrostatic field applied perpendicular to 
the nanotube axis [23,[2l|. In both cases the exciton prop- 
erties are mediated by collective plasmon excitations in 
CNs [13, [2^. In the latter case (Fig. [T|), we have shown 
recently that the QCSE allows one to control the exciton- 
interband-plasmon coupling in individual undoped CNs 
and their optical absorption properties, both linear [2^ 
and nonlinear [2l|, accordingly. 

In general, plasmons cannot be excited by light in op- 
tical absorption since they are longitudinal excitations 
while photons are transverse. In small-diameter (~ 1 nm) 
semiconducting CNs, light polarized along the CN axis 
excites excitons which, in turn, can couple to the near- 
est (same-band) interband plasmons [l0lllll.[20j. Both of 
these collective excitations originate from the same elec- 
tronic transitions and, therefore, occur at the same en- 
ergies ^ 1 eV, as opposed to bulk semiconductors where 
they are separated by tens of eVs. They do have different 
physical nature. Their coexistence at the same energies in 
CNs is a unique general feature of confined quasi- ID sys- 
tems where the transverse electronic motion is quantized 
to form ID bands and the longitudinal one is continuous. 

The formation of coupled exciton-plasmon excitations 



can be viewed as an additional non-radiative channel (in 
addition to phonons (l^l and defects [l^) for the exciton 
relaxation in CNs, where optically excited excitons decay 
into low-energy interband plasmons. In so doing, excitons 
generate the quanta of plasma oscillations on the CN sur- 
face, on the one hand, and this shortens their lifetime, 
on the other. Thus, by varying the exciton-plasmon cou- 
pling strength using the QCSE one controls both the ra- 
diative emission from an individual CN and surface elec- 
tric field fluctuations associated with plasmons generated 
by excitons on the CN surface. This latter phenomenon 
is pretty much similar to the SPASER effect (Surface 
Plasmon Amplification by Stimulated Emission of Radi- 
ation) reported earlier for hybrid metal-semiconductor- 
dielectric nanosystems [2^, and is the focus of the stud- 
ies here for individual small-diameter CNs. The non- 
radiative exciton-to-plasmon energy transfer is shown to 
result in the large population numbers of coherent surface 
plasmons associated with high-intensity coherent oscillat- 
ing fields concentrated at nanoscale across the CN diam- 
eter along the CN length. These strong local fields can 
be used in a variety of new CN based optoelectronic ap- 
plications, including near-field nonlinear-optical probing 
and sensing, optical switching, enhanced electromagnetic 
absorption, and materials nanoscale modification. 
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(in part from Ref. [20| , wherein similar data for other CNs 
are found). In (b), we see the anti-crossing ~ 0.1 eV, both 
in the energy-momentum plane and in the energy-field 
plane, revealing the QCSE as an efficient tool to control 
the exciton-plasmon coupling strength in individual CNs. 

When the exciton is excited and the CN surface 
EM field subsystem is in the vacuum state, the time- 
dependent wave function of the whole system "exciton -|- 
surface EM field" is of the form 

= E^/(k>Oe-'^^('')*/''|{l/(k)})e.|{0}) (1) 

+ ^ du;C{k,u;,t)e-'-*\{0}),,\{l{k,u;)}). 
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FIG. 2: (Color online) Calculations for the (11,0) CN. 
(a) Fragment of the energy dependence of the dimensionless 
(normalized by /2nfi) axial surface conductivity azz along 
the CN axis [2y|. Ovals mark exciton {En, E22) and inter- 
band plasmon (Pii, P22) excitations [peaks of Yiea^z and 
Yie(l/(Jzz), respectively], (b) Exciton-plasmon dispersion re- 
lation as a function of the perpendicular electrostatic field ap- 
plied (see Fig. [T]) and longitudinal momentum for the lowest 
bright exciton [En in (a)] when coupled to the nearest inter- 
band plasmon [Pn in (a)]. Dimensionless momentum and di- 
mensionless energy are defined as 'ib[M omentum] / {2iTh) and 
[Energy] /2'yo, respectively, where b — 1.42 A is the C-C inter- 
atomic distance and 70 — 2.7 eV is the C-C overlap integral. 



II. PLASMON GENERATION BY EXCITONS 

In small-diameter semiconducting carbon nanotubes, 
because of their quasi-one-dimensionality, excitons are 
excited by the external electromagnetic (EM) radiation 
polarized along the CN axis [27j. As a consequence, they 
have their transition dipole moment and translational 
quasi-momentum both directed along the nanotube axis 
{z axis of the problem — see Fig. [T] cylindrical coor- 
dinates are used). That is why they are able to cou- 
ple to their neighboring longitudinal interband plasmon 
modes [ll|, HO] • Figure [5] shows an example for the exci- 
tons and interband plasmons [Figj^Ka)], and the exciton- 
plasmon dispersion relation [Figl^^b)] in the (11,0) CN 



Here, |{l/(k)})ea; is the single-quantum Fock state with 
one /-internal-state exciton excited of quasi-momentum 
k — {k^,kz} with quantized component (transverse 
quantization that results in electron-hole subbands) and 
continuous kz component representing the longitudi- 
nal motion of the exciton. The excited single-quantum 
state |{l(k, cij)}) represents longitudinally polarized sur- 
face EM mode of frequency w (the plasmon). The re- 
spective vacuum states are |{0})ea; for the exciton sub- 
system and |{0}) for the surface EM field subsystem. 
The coefficients C/(k,t) and C(k, lli,^) stand for the 
population probability amplitudes of the exciton sub- 
system and field subsystem, respectively. The exciton 
energy is of the form i?/(k) ~ Ef{'k) — ih/r, where 

Ef{k) = Eiil{k^) + h^kl/2M^^{k^) stands for the total 
energy of the /-internal-state exciton and r represents 
the phenomenological relaxation time to account for all 
possible (but into plasmons) slow exciton relaxation pro- 
cesses (normally attributed to the exciton-phonon scat- 
tering [2^ lU)- The energy i?/(k) consists of the ex- 



citon excitation energy Eill{k^) 



Eg(k^) + Eif\k^), 



where Eg{k^p) — Se^k^p) + eh{kp) is the band gap with 
£e,h being the transversely quantized azimuthal electron- 
hole subbands, e\J'^ is the (negative) binding energy of 
the exciton, and the kinetic energy of the translational 
longitudinal motion of the exciton with the effective mass 
Mex = iTie + irih, where rUe and ruh are the (subband- 
dependent) electron and hole effective masses. 

The time dependent Schrddinger equation with the 
Hamiltonian of the whole system " exciton -I- surface EM 
field" applied on the wave function ([Ij results in the 
set of the coupled differential equations for C/(k, i) and 
C(k,w,t) as follows [20] 

ihCf (k, t) e-*^/(k)*/?' ^ duj g^+^ (k, k, uj) C(k, UJ, t)e-'^\ 
ih C7(k, c., t ) * = ^ [g}+) (k, k, c.)] * Cf (k, t)e- -^^f C^)*/''. 

Here, the interaction matrix element squared is given by 



gi^+^(k,k,c.) 



;i^o]rf{]'c^3 



Re 



1 



o-zz(k,a;) ' 



(2) 
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where 4 = En(0|(dn)z|/) with (0|(d„),|/) being the 
electric dipole transition matrix element where the /- 
internal-state exciton is excited in the lattice site n of the 
CN surface, azz stands for the CN surface axial conduc- 
tivity with Ile{l/azz) representing the plasmon density 
of states [DOS, see Fig. [2] (a)], Rcn is the CN radius, 
and 5*0 = (3y3/4)&2 is the area of an elementary equi- 
lateral triangle selected around each carbon atom in a 
way to cover the entire CN surface, b= 1.42 A is the C-C 
interatomic distance. 

In terms of the probability amplitudes Cf(k,t) and 
C(k, the exciton emission intensity distribution to 
generate plasmons is given by the final state probability 
at long times corresponding to the decay of all initially 
excited excitons, i.e., /(k,cj) — |C(k,w,i— >-(X))p. This is 
related to the exciton probability amplitude C/(k, t), ac- 
cording to the equations above, and thereby is also asso- 
ciated with light absorption by excitons. The peak inten- 
sity represents the long-time plasmon population gener- 
ated by optically excited excitons with the momentum k. 



III. PLASMON INDUCED SURFACE FIELD 

Plasma oscillations generated by the non-radiative ex- 
citon decay into the nearest interband plasmon mode 
[Fig. [5] (a)] can be viewed as standing charge density 
waves due to the periodic opposite-phase displacements 
of the electron shells with respect to the ion cores in the 
neighboring elementary cells on the CN surface. Such 
periodic displacements induce local coherent oscillating 
electric fields of zero mean, but non-zero mean-square 
magnitude, concentrated at the nanoscale across the 
diameter throughout the length of the nanotube. The 
mean-square longitudinal local field magnitude can be 
calculated as the observable, the expectation value 



i?2(n)-([-V„^(n)]^ 



(3) 



of the quantum electrodynamical (QED) longitudinal- 
electric-field operator at the lattice site n = {i?cw, z„} 
of the CN surface in the state with the wave function ([IJ . 
This results in (see Appendix A) 



fiL 



(4) 



k — kin , k ^ 



f°° 1 

X / dujuj^Re 

7o azziK^) 



iV(k,^)+- 



where 7V(k, w) — (47rc/L)/(k, a;) is the number of plas- 
mons generated by excitons with the momentum k of the 
1st Brillouin zone on the surface of the tubule of length L, 
and fexiKT) = exp[-h'^Pj2Me4k^)kBT]/Qex is the 
exciton momentum distribution function with the parti- 
tion function QexiT) = J2kexp[-?fPj2Mex{k^)kBT]. 

Equation (|H) tells that in order to produce strong lo- 
cal mean-square fields, the exciton non-radiative emis- 
sion resonance [given by /(k, lo Eexc)] must overlap 



(11,0) 
T=10 K 




FIG. 3: (Color online) Low (a) and high (b) longitudinal- 
momentum-averaged plasmon population (arb. units) for the 
non-radiative decay of the 1st bright exciton in the (11,0) CN. 
See Fig. [5] caption for the dimensionless energy. 



with the energy Ep of the neighboring plasmon DOS res- 
onance {the peak of Re[l/(Tzz{i^ ^ Ep)]}. By way of ex- 
ample of the (11,0) CN in Fig. [51 we see that this can be 
achieved by carefully tuning the exciton excitation en- 
ergy Eexc and the nearest interband plasmon energy Ep 
by means of the QCSE. As this takes place, both E^xc 
and Ep experience red shift with increasing perpendicular 
electrostatic field. However, the red shift of E^xc is very 
small [barely seen in the constant momenta energy-field 
planes in Fig. ^ (b)] due to the negative field- dependent 
contribution of the exciton binding energy Eb [20| • As a 
consequence, E^xc and Ep [En and Pu in Fig. [2] (a)] ap- 
proach as the field increases, pushing the coupled exciton- 
plasmon system into the strong coupling regime where all 
of the optically excited excitons decay non-radiatively to 
generate interband (same-band) plasmons, yielding the 
peak optical absorption at the same time. As the temper- 
ature grows up, higher momenta excitons start contribut- 
ing to the process, lowering the field necessary to achieve 
the strong coupling regime [blue contrast in Fig. [5](b)]. 

Figure [3] shows plasmon population numbers aver- 
aged over the longitudinal momentum distribution f28j . 
Sfc-/ea:(k, T)iV(k, w), also representing light absorption 
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FIG. 4: (Color online) Local surface field amplitude, E^, as 
a function of temperature and perpendicular electrostatic field 
applied, given by Eq. Q for plasmons generated by the non- 
radiative decay of the 1st bright exciton in the (11,0) CN. 



by excitons, calculated at low and high temperatures for 
the first bright exciton in the (11,0) CN with E^.^^ and Ep 
[E'li, Pii in Fig. [5] (a)] tuned using the QCSE. The same 
exciton-plasmon parameters and their perpendicular- 
electrostatic-field dependences are used as in Ref. j23| 
{Eexc = 1.21 eV, Eb = -0.76 eV, Ep = 1.51 eV in the 
zero field, and r = 30 fs to account for the exciton- 
phonon relaxation [24|). We see the dramatic increase in 
the peak intensities, associated with increased optical ab- 
sorption, when the field strength exceeds 5 x 10^ V/m, 
both at low and at high temperatures [1^. Rabi split- 
ting occurs as the field drives the exciton-plasmon sys- 
tem into the strong coupling regime, whereby the effec- 
tive plasmon generation starts. Temperature generally 
smoothens the effect due to higher momenta excitons 
contributing to the process. 

Large plasmon population numbers are, according to 
Eq. Q, associated with coherent oscillating fields con- 
centrated locally throughout the CN surface. Figure |4] 
shows the calculations for the (11,0) CN. Huge local sur- 
face fields ~ 10^ V/m, just a few orders of magnitude 
less than intra-atomic fields, are created under the reso- 
nance conditions. The effect slightly decreases with tem- 
perature, but it starts at lower electrostatic fields due 
to higher momenta excitons contributing to the plasmon 
generation. Strong local surface fields created are the re- 
sult of the efficient energy conversion, whereby the ex- 
ternal EM radiation energy absorbed to excite excitons 
transfers into the energy of high-intensity coherent local- 
ized optical-frequency fields of charge plasma oscillations. 



IV. CONCLUSIONS 

The effect presented here for individual single wall CNs 
is analogous to the SPASER effect reported earlier for hy- 



brid metal-semiconductor-dielectric nanostructures [26|. 
The effect is universal as it originates from the trans- 
verse quantization of electronic degrees of freedom in 
quasi-lD systems. It can manifest itself in densely packed 
aligned nanotube films as well, both through plasmon 
enhanced inter-tube Casimir interactions [N = 0-term 
in Eq. Q], as it is recently demonstrated for double wall 
CN systems [31J, and through the exciton-to-plasmon en- 
ergy transfer tuned by means of the QCSE. In the latter 
case, plasmon-induced coherent local surface fields can 
be used in a variety of new tunable optoelectronic ap- 
plications with both individual nanotubes and nanotube 
composites, such as enhanced electromagnetic absorption 
and optical switching, near-field nonlinear-optical prob- 
ing and sensing, materials nanoscale modification. 
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Appendix A: Derivation of Eq. ((4]) 

Here, the exciton-plasmon coupling theory of Ref. j20| . 
developed earlier for individual single wall semiconduct- 
ing CNs, is used to derive Eq. Q for the mean-square 
longitudinal local field induced by plasmons generated by 
optically excited excitons on the CN surface. In this the- 
ory, the nanotube is modelled by the neutral, infinitely 
long, infinitely thin, anisotropically conducting cylinder, 
where only the axial conductivity a^z is taken into ac- 
count while the azimuthal one is neglected being strongly 
suppressed by the transverse depolarization effect. The 
vector ez of the orthonormal cylindrical basis {e,-, e^p,ez} 
is directed along the CN axis as shown in Fig. [TJ 

The longitudinal surface electric field operator at the 
lattice site n = {Rcn, fmZn] of the CN surface is of the 
form (Gaussian units) 

f°° - II 

- Vn<^(n) =/ dwE (n,w) + /i.e., (Al) 
Jo 

where 

E"(n,a;) II G(n,m,tj) -5(111,0;) (A2) 

m 

is the longitudinally polarized Fourier-domain surface 
electric field operator in the Schrodinger picture with 
llG(n, m,a;) being the longitudinal part of the classical 
EM field Green tensor on the CN surface, and 

T/ N /^wRecrzz(w) - 

J(n,cj) = y /(n,w)ez (A3) 

representing the CN axial surface current density opera- 
tor defined in such a way as to preserve the fundamental 
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QED equal-time commutation relations for the EM field 
operator components in the presence of the CN medium 
assisted absorption (see, e.g., Ref. [13]). Here, f{n,uj) 
along with its counter-part /''' (n, w) stand for the scalar 
bosonic field operators to annihilate and create, respec- 
tively, single-quantum EM field excitations of the fre- 
quency Lu at the lattice site n of the CN surface. 

The CN axial current density operator convoluted with 
the longitudinal EM field Green tensor per Eq. (jA2|) 
yields the correctly quantized longitudinal surface elec- 
tric field operator (jAll) . whose z-component associated 
with the zz-component of the longitudinal Green tensor 

\\r (r, m ,A- Cy^d{ipn-ip„,) 

"i^zz[n,m,uj) — — ^ ' 

SnRcN a-zz[io) 

(A4) 

is responsible for plasmon-induced coherent local surface 
fields as given by Eqs. ([3]) and (|4|). 

The longitudinal Green tensor component (|A4|) has the 
following property 



"GL(n>m:^)"Gz2(n, m',tj') = 

n 

m 



a ^m' ) / ^ 



(jj— Ul' 



32t:'^ RcNcr*^ {uj)a^^ (w') 



(A5) 



This can be proven replacing the summation over lat- 
tice sites by the integration over the nanotube surface 
according to the rule 



^ So 



1 

— / dipriRcN 
JO Jo 



— oo 



where So = (3%/3/4)62 (6= 1.42 A being the C-C distance) 
is the equilateral triangle area selected around each car- 
bon atom in such a way as to cover the entire CN surface. 
The integration over Lpn is trivial, and the one over z„ is 
performed by representing the cosine-function in the ex- 
ponential form, writing the integral as follows 



dzn-.. = lim 



L/2 
-L/2 



with L being the CN length, and dividing it into parts 
by means of the equation 



^iu:\Zn-Zrr,\/c 



iuj{zji — Zm) /c 



'iuj{Zn~ Zm) I C 



This, after some algebra with the use of the following 
representation for the ^- function (see, e.g., Ref. (ssj ) 



(5(x) — — lim 

t — > +00 



(A6) 



brings one to the equation 

r+L/2 



lim 

L— > + CXD 



-L/2 



= 2c cos 



(^m ^m') 



2ttS{u;- uj')- 



to eventually result in Eq. (jASp . The (5-function repre- 
sentation (|A6|) can also be written in a more general form 



Sj.)^'"!^:'}- hm ^" 



p= 1,2,3,... (A7) 



2ntP-^ ^^^0 {x + ie)P 

which is proved by using the complex plane contour in- 
tegration technique and the Cauchy integral theorem. 

The property (jA5[) can be converted to the momentum 
space to yield 



^llG:,(n,k,w) llG,,(n,k',L.') 



Nc^ 6uu' 



6ATr^RcNat,iuj)azziLo') 



(A8) 



where 



II G.. (n, k, w) - ^ ^ II G.. (n, m, to) e^ "^ 

with N = 2ttRcnL/ So being the number of the lattice 
sites (carbon atoms) on the CN surface, and 

1 Z"^'' 



X lim — 

L— f-t-oo L 



-L/2 



L/2 



(A9) 



From Eqs. (jAll) - (jA4[) . taking only the single-quantum 
states into account [in line with Eq. ([T])], one has 



poo 

-V„(^(n)]^= dudu'A*{u)A{uj') 
Jo 

^llGL(n,k,c.) llG,,(n,k>') 

k,k' 

2p{k,uj)f{-k',Lj') + Jk,k'<5(w - uj') 



(AlO) 



where 



A{uj) = I — Y ' 



^A*(w), (All) 



/t(k,a.) = ^5]/t(m,a.) 

^ m 

/(k,w)=[/t(k,,,) 



with the k summation running over the first Brillouin 
zone of the CN. Further, in view of the fact that all the 
lattice sites are equivalent, one has 



hV„^(n)]2 = l^hV^^(m)] 



(A12) 
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for any arbitrary n. Equations (IA12[) and (jA8l) allow one 
to reduce Eq. (jAlOp to the form 

[-V„^(n)]^= dujduj'A*{uj)A{uj') (A13) 

^0 



32tt^ RcNcrl^ (a;)(T^z(iu') 



2TrS{uj-uj')- 



E 



Now, the expectation value ([3]) can be calculated. 
Equation ^ contains two types of averaging. They are 
the quantum mechanical averaging with the wave func- 
tion ([T} and the statistical one over the exciton quasi- 
momentum k. Since the typical CN exciton excitation 
energies are > 1 eV (see, e.g., Ref. [l^), the excitons 
are assumed to be originally excited by the external EM 
radiation. Once excited, the thermal exciton momentum 
relaxation starts in subbands. This situation is properly 
accounted for by the use of the density operator 



(A14) 



where k = {fci^, fc^} runs over the first Brillouin zone of 
the nanotube, 

|V'(k,t))^^Q(k,i)e-^^(*'^*/''|{l/(k)})e.|{0}) 
/ 



duC{k,w,t)e-'"'\{0}),,\{l{k,uj)}) (A15) 



with the normalizing condition (?/;(k, i)|?/j(k', t)) = '^kk' 
[cf. Eq. (P)], and 



fexiKT) = 75— exp 



21,2 



(A16) 



is the normalized classical (Boltzmann) exciton momen- 
tum distribution function with the partition function 



Qex{T) = y^exp 



2Mex{k^)kBT 



(A17) 



This latter one can be evaluated explicitly as follows 



Qex{T) = 7;- / dk^ exp 



21,2 



2Mex{k^)kBT 



ksT 
2T:h'^ 



^ JMe,(fc^)Erf 



2i,2 



2Mexik^)kBT 



mL 



M,,.kuT 



21:%'' 



■Erf 



h'^kl 



2MexkBT 



Here, the first Brillouin zone of the CN of (m, n) type 
{n < m) is taken to be consisting of m parallel lines, as 
per quantized kip = fci^(s) — s/Rcn with s — l,2,...,m 
and RcN = {\/3b/ 271)^/111^ + mn + rfi , each of length 
2k^ = 2B/kip{m), where 2B = 2{Att'^ /i^/Zb'^) is the rect- 
angular area of the reciprocal space covered by the lines. 
This yields k^ = 27r/3& for the (m,0) type CNs (zigzag) 
and kz = 27r/-\/3& for the {m,m) type CNs (armchair), 
in particular. The last expression above is the approx- 
imation neglecting the difference in the effective masses 
for different subband excitons, yielding M^^ikip) w M^x 
for all A:^, to result in X]fe,^{---} ~ {•••IS^Lil = {...}m. 
Using the density operator (IA14[) . one has 



£;f(n)= lim Tr{[-V„^(n)]2p(t)}. 



(AI8) 



This expands further as follows 

i?2(n) = lim ^(^(k,t)|[-V„^(n)]2p(i)|V;(k,0) 

k 

hm /- T) (V(k, t) I [- Vn (^(n)]2 |^(k, 0) , 



wherein, upon substitution of Eqs. (jA13[) and (jA15[) in 

it, one has 

(^(k,i)lE/^(k'.^)/(k''^')l^(k,0> 



C*(k,w,t)C(k,w',i)e 



— i{uj' —uj)t 



to result in 



El{n)= Van fex{KT) duj duj' A* {uj)A{uj') 



32'K^RcN o-lzi^)cTzziuj') 



(A19) 



2nS{uj-uj')- 



C* (k, cj, t)C(k, Lo', t) e-*('^'"'^)* + ^6{ijj'- J) 



Multiplying the square brackets in here yields four terms 
to be taken in the limit of < — > +00, followed by the 
integration over lj' . The four terms are 

2^C* (k, cj, t)C{k, J, t) e~'("'-")*(5(cj- J) 



-iC*(k,w,t)C(k,w',i)- 



+ ttS{uj— uj')5{uj— lj') 

If/ /N * 

+ 5{u:~uj')- . 

2 cj' — cj 

The first term here is trivial to deal with. The second 
term is treated by formally adding a positive imaginary 
infinitesimal to its denominator to write uj' — uj -\- it with 
e— )--|-0, followed by using Eq. (IA6|) . The term sums up 
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with the first one. To treat the third term correctly, one 
first has to represent one of the (5-functions as follows 

1 /•+i/2 

S{uj - uj') = lim / dz„ e-'("-" (A20) 

L^+oo 27rc J_L/2 

and then integrate the whole contribution over w'. This 
gives the term proportional to the CN length L. Finally, 
the last, fourth term can be shown to yield zero principal 
value by noticing that it is an odd function of its variable. 
The (5-function makes it zero for all w ^ w', where it must 
be zero, too, since it is odd. More accurate proof of this 
fact can be done in several ways, e.g., by using the iden- 



tity d'{x) = —S{x)/x with x = u) — uj', followed by partial 
integration of the whole contribution, or starting from 
the (5- function representation (jA20p . followed by manip- 
ulating the integration order over uj and cj' with the use 
of the general properties of the integrands. 

Inserting the contributions discussed into Eq. (|A19p . 
adding them all together in the limit of t — > +oo with 
the use of Eq. (jAlip for the coefficient A{uj), and taking 
into account the fact that Re[l/crz2(a;)] is a sharp peak 
structure [where peaks should be attributed to different 
kip{s) = s/RcN, s = 1,2,..., to represent different inter- 
band transitions [l^; see an example in Fig. [5] (a)], one 
eventually arrives at the equation (|4]). 
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